MODELING THE DISRUPTION OF THE GLOBULAR CLUSTER PAL 5 BY GALACTIC TIDES 



^1" 
O 
O 

(N 

C 



(N 



> 
(N 
(N 

o 
o 

Or 

o ■ 



Walter Dehnen 123 , Michael Odenkirchen 3 , Eva K. Grebel 43 and Hans-Walter Rix 3 

walter.dehnen@astro.le.ac.uk; odenk@mpia.de; grebel@astro.unibas.ch; rix@mpia.de 

ABSTRACT 

The globular cluster Pal 5 is remarkable not only because of its extended massive tidal tails, but also for its very low 
mass and velocity dispersion and its size, which is much larger than its theoretical tidal radius. In order to understand 
these extreme properties, we performed more than 1000 A-body simulations of clusters traversing the Milky Way on the 
orbit of Pal 5. Tidal shocks at disk crossings near perigalacticon dominate the evolution of extended low-concentration 
clusters, resulting in massive tidal tails and often in a quick destruction of the cluster. The overlarge size of Pal 5 can 
be explained as the result of an expansion following the heating induced by the last strong disk shock ~ 150 Myr ago. 
Some of the models can reproduce the low observed velocity dispersion and the relative fractions of stars in the tails and 
between the inner and outer parts of the tails. Our simulations illustrate to which extent the observable tidal tails trace out 
the orbit of the parent object. The tidal tails of Pal 5 show substantial structure not seen in our simulations. We argue that 
this structure is probably caused by Galactic substructure, such as giant molecular clouds, spiral arms, and dark-matter 
clumps, which was ignored in our modeling. 

Clusters initially larger than their theoretical tidal limit remain so, because, after being shocked, they settle into a 
new equilibrium near apogalacticon, where they are unaffected by the perigalactic tidal field. This implies that, contrary 
to previous wisdom, globular clusters on eccentric orbits may well remain super-tidally limited and hence vulnerable 
to strong disk shocks, which dominate their evolution until destruction. Our simulations unambiguously predict the 
destruction of Pal 5 at its next disk crossing in ~ 1 10 Myr. This corresponds to only 1 % of the cluster lifetime, suggesting 
that many more similar systems could once have populated the inner parts of the Milky Way, but have been transformed 
into debris streams by the Galactic tidal field. 
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1. INTRODUCTION 

Amongst the globular clusters of the Milky Way galaxy, Pal 5 is 
exceptional in many respects. First, it is one of the faintest and 
least massive of these objects; its velocity dispersion is so small 
that it could only recently be determined, using high-resolution 
spectroscopy, to be below lkms -1 (Odenkirc hen et aii r2002, 
hereafter paper I). Second and most remarkable, it has been de- 
tected, using d ata from the Sloan Digital Sky Survey (SDSS; 
I York et all 120001) . to possess a pair of tidal tails, which ex- 
tend at least 4 kpc in opposite directions fr om the cluster and 
conta i n more stars then the cluster it self ( Odenkirchen et alJ 
1200 It iRockosi et alJl2002t lOdenkirchen et alJl2003l hereafter 
paper II). Pal 5 is actually the first and only globular cluster 
so far for which such tails have been detected at a compa- 
rable level of significance. Thus, this is the first direct ev- 
idence that mass-loss induced by the Galactic tidal field can 
substantially affect the evolution of globular clusters, a process 
which is speculated to have destroyed many low-mass clusters 
of an initially much rich er Galactic system of clusters (e.g., 
IVesperini & Heggiell997l) . 

The tidal field generated by a smooth Galactic mass distri- 
bution stretches the cluster in a direction towards the Galactic 
center and weakly compresses it in tangential directions. The 



stretching creates drains through which stars from the outer 
parts of the cluster are carried away. This limits the bound 
part of a cluster which is in equilibrium with the tidal field to a 
tidal radius, r t ;d, whose size may be estimated by equating the 
outwardly directed stretching tidal force to the cluster's inter- 
nal gravitational attraction. For a cluster of mass M orbiting at 
galactocentric radius R in a galaxy with constant circular speed 
v c , we get 5 
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For eccentric orbits, r tl d varies along the orbit, and is smallest 
at perigalacticon. The strong tidal force at perigalacticon acts 
only for a short fraction of the orbital period, and passages of 
the perigalacticon are often referred to as tidal shocks or 'bulge 
shocks'. 

The tidal field generated by the stellar disk with its steep ver- 
tical gradient is different. It is much stronger than that of the 
smooth halo and bulge, and it is compressive. The latter is be- 
cause, whilst the cluster is crossing the disk, its stars feel the 
additional attraction of the disk stars within the cluster. How- 
ever, the strong compression only acts for a short time. This 
allows to use the impulse approximation, which gives for the 
velocity chan ge of a star at vertical distance z from the cluster 
center ( Ostriker, Spitzer & Chevalierll972t ISpitzeJ 19871) 
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where V is the vertical velocity of the cluster crossing the disk 
and g m w 2-7T G T,(R) the maximal vertical acceleration exerted 



5 In the literature this equation often comes with an extra factor 1/2, which 
occurs either if the orbit is assumed to be circular or the Galactic potential to be 
that of a point mass. If both assumptions are made, the factor is 1/3. 
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by a (thin) disk with surface density These changes in 

velocity result in relative changes of the specific energies of the 
stars of the order of 

(jAgj) crg m rh 2ff m r h 
(\E\) ~V<|£|)~ ff V ^ 

where we have used sw rh/2 with r*h denoting the cluster 
half-mass radius. Low-concentration systems tend to have large 
half-mass radii and low velocity dispersions, which makes them 
very susceptible to disk shocks, in particular at small galacto- 
centric radii where and hence g m are large. Not surpris- 

ingly therefore, low-concentration clusters are hardly found in 
the inner Galaxy, where not only the shocks are much stronger 
but also occur much more frequently than for orbits in the 
Galactic halo. 

The velocity changes also res ult in an average heating per 
star o f (AE) « 2gl l r 2 i /W 2 (Ostriker. Snit zer & Chev alier 
119721 ISpitzerlfl987h . One may use this result to estimate the 
"shock heating time" t 8 h, the time scale over which disk shocks 
significantly affect the cluster dynamics: 

-(E) _ 3 a 2 V 2 
tsh ~ tdisk (AE)~ tdisk Aglrl C) 

(iGnedin. Lee & Ostriker 1999). Here, idisk is the time between 
strong disk shocks, which for eccentric orbits, such as that of 
Pal 5, equals the orbital period P. 

Apart from tidal forcing, the evolution of globular clusters 
is also driven by internal processes, such as stellar mass loss, 
two-body relaxation and binary interactions. These internal pro- 
cesses cause the cluster to (eventually) undergo a core collapse 
and result in mass-segregation and evaporation of preferentially 
low-mass stars. While the evolution of isolated globular clus- 
ters has been the subject of many studies, the combined effects 
of these int ernal and the external proc esses have been rarely in- 
vestigated. IGnedin & Ostrikei] ill 9971) used Fokker-Planck sim- 
ulations including two-body relaxation, tidal limitation and disk 
and bulge shocks to investigate the dissolution of the individ- 
ual Galactic globular clusters (though their adopted mass for 
Pal 5 is five times larger tha n our best curre nt value). With the 
same method, Gnedin, Lee & Ostrikerl l[l999) studied the evolu- 
tion of globular clusters with various concentrations c, defined 
as c e log(ru m /r core ) from the cluster's limiting and core ra- 
dius, and ratios (3 = £ r hAsh between the two-body relaxation 
time i r h and the shock heating time (|3}. They considered clus- 
ters with c e [0.6, 2.6] and (3 <E [1CT 5 , 10 2 ] orbiting the Milky 
Way on the orbit of the cluster NGC 6254 and concluded that (i) 
for j3 > 0.1 disk shocks dominate the cluster evolution and may 
lead to quick destruction, (ii) for smaller values of (3 tidal forc- 
ing accelerates the two-body-relaxation driven evolution, but 
also that (iii), in response to disk shocking, the cluster com- 
pacts, substantially reducing (3 and diminishing the further im- 
portance of disk shocks. Pal 5 has c ~ 0.6 (paper II) and (3 ~ 10 
("section l2~4l . which means that, according to Gnedin et al., the 
evolution of Pal 5 is entirely dominated by disk shocks. 

Collisional A-body-simulation studies which investigate in- 
ternal processes alo ng with a time varying tidal field either 
ignore disk shocks (Baumgardt & Makino 2003) or account 
for disk shocks but assum e an otherwise constant tidal field 
dVesperini & He ggiel 19971) . 

Unfortunately, neither of these studies is directly applicable 
to Pal 5, since they either exclude relevant external processes or 



do not cover the globular-cluster parameters relevant for Pal 5. 
Gnedin et al., for example, assumed in their study that disk 
shocks are slow and rare, i.e. the cluster-internal dynamical time 
is shorter than the duration of the shock and much shorter than 
idisk, whereas for Pal 5 shocks are fast and frequent (see sec- 
tion l2.4> . Also, the orbit of Pal 5 differs significantly from that 
of NGC 6254, as adopted by these authors. 

Moreover as far as we are aware, all studies of globular clus- 
ter evolution start with clusters limited by their (perigalactic) 
tidal radius, assuming that this limitation is quickly achieved by 
the tidal force field. As we will see below, this assumption is not 
justified, neither theoretically as our simulations will demon- 
strate, nor observationally, since Pal 5 is an excellent counterex- 
ample, see section l2~3l 

The primary goal of the present paper is to study, via detailed 
A-body simulations, the effect of Galactic tides on a globu- 
lar cluster moving on the orbit of Pal 5 and to quantitatively 
compare the observable properties of the models with those of 
the cluster in an attempt to understand its current dynamical 
state and to constrain the formation history of this object. As 
a byproduct, our simulations will also give valuable insights 
into the evolution of low-concentration globular clusters with 
(3^1 experiencing strong disk shocks in conjunction with a 
time varying tidal field. 

In section 13 we summarize the observed structural and dy- 
namical properties of Pal 5. The A-body simulations are pre- 
sented in section The dynamics of the simulated tidal tails 
are discussed in detail in section |4] In section |5] we compare 
our simulations directly to the data for Pal 5. The results of the 
A-body simulations and the implication for Pal 5 are discussed, 
respectively, in sections[6]and[7]and summarized in section|8] 

2. THE OBSERVED PROPERTIES OF PAL 5 

We shall now summarize the structural and kinematic properties 
of the globular cluster Pal 5, which are relevant for the present 
study. Most of these properties have been presented and are 
discussed in papers I & II. 

2.1. The Orbit 
As shown in paper II, the orbit of Pal 5 is rather tightly con- 
strained by its position and radial velocity in conjunction with 
the observed orientation and curvature of the tidal tail. This is 
the case, because the tail stars deviate from the cluster orbit only 
very little, as is evident, for instance, for the simulated tidal tail 
in Figure |6] below. Figure [2 shows the meridional projection 
of that orbit obtained with our standard model for the Galac- 
tic gravitational potential (see section l3~TI for details). Differ- 
ent assumptions for the Galactic potential, e.g. different circular 
velocities, result in very similar orbits, except for their periods 
(paper II). In particular, the perigalactic and apogalactic radii 
are always at about 5.5 and 19kpc. 

2.2. Mass and Velocity Dispersion 
With a total luminosity of only M v = -4.77 ± 0.20 (paper I) 
Pal 5 is one of the faintest Galactic globular clusters. It also has 
an un usually flat stellar luminosity function (Grillmair & Smith 
1200 ll) implying that its mass-to-light ratio is atypically low. 
Together, this let us estimate in paper I the total mass to be 
5200 ± 700 Mq, substantially less than previous estimates. 

In paper I, we measured the line-of-sight velocity dispersion 
trios from high-resolution spectra of 18 giant stars within 6' of 
the cluster center to be at most 1.1 km s _1 . However, the line- 
of-sight velocity distribution is significantly non-normal but has 
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Fig. 1. — Orbit of Pal 5, projected into the meridional plane, for the assumed 
Galactic potential. The current position of Pal 5 is indicated by the star; the 
trajectory is plotted for the last 2.95 Gyr (solid), corresponding to 10 radial pe- 
riods, and the future 0.5 Myr (dotted). 

extended high-velocity tails, which dominate the calculation of 
(j\ os . When modeling this distribution as Gaussian for the clus- 
ter dynamics plus a contribution from binaries to account for the 
tails, we obtained for the dynamically relevant velocity disper- 
sion er los = 0.22;°;}o km s _1 . 

As we discussed in paper I, these estimates for the mass and 
velocity dispersion are consistent with the best-fit King model 
and the assumption of dynamical equilibrium. However, this 
does not imply that the cluster is in equilibrium and we shall 
now see that it cannot possibly be. 

2.3. Extent and Tidal Radius 
Figure[2]shows the surface density, averaged over annuli, of the 
cluster and the tidal tail. The profile of the cluster itself (which 
can easily be obtained from data in directions perpendicular to 
the tail) appears to be truncated at a limiting radius of about 
16'. 1 « 107 pc, see Fig. 5 of paper II. This radius must be com- 
pared to the theoretical tidal radius (equation Q for the cluster, 
which for v c — 220 km s -1 and the mass estimate given above 
is ~ 54 pc at the current location of Pal 5 and only ~ 24 pc at the 
perigalacticon of its orbit. This latter value is much smaller than 
run, and actually equals the cluster's present-day core radius. 

Thus, the cluster is much larger than its tidal radius and 
cannot possibly be in equilibrium with the Galactic tidal force 
field 6 . This is surprising and shows that the current dynamical 
state of Pal 5 must be a peculiar one. In particular, any estima- 
tion which is based on the assumption of equilibrium may well 
be in error. 

2.4. The Importance of Disk Shocks 
After having assessed the various properties of Pal 5 and its or- 
bit, we may now consider the importance of disk shocks at disk 
crossings. The tidal force acting on a star at position r with 



6 In previous studies, this discrepancy was not detected because (i) the mass 
of the cluster was overestimated by a factor of 4-6 and (ii) the perigalactic radius 
was unknown. 
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Fig. 2. — Radial profile of the azimuthally averaged surface number density of 
SDSS stars in Pal 5 and its tails as measured in paper II (top), and cumulative 
number distribution averaged over both tails (bottom). The typical mass of Pal 5 
stars in SDSS is about 0.8 Mq. 

respect to the cluster center is given by 

F tid (r,t) = (r-V)V$(R(t),t) (4) 

where $ is the gravitational potential of the Milky Way and 
R(t) the galactocentric position of the cluster center at time t. 
Equating LFtidl to the cluster internal attraction GM/r 2 gives 
relation Q for the radius r if $ is taken to be the potential of a 
singular isothermal sphere. 

In Figure|3j we plot the strength of the tidal force field, quan- 
tified by the eigenvalues of 9 2 $ / 'dx. L dxj as function of time on 
Pal 5's orbit in our model for the Galactic potential (see section 
13. 11 . On top of a smooth underlying tidal field, whose dom- 
inant effect is a stretching (the absolute largest eigenvalue is 
negative, i.e. dotted in Fig. |3jl and whose strength varies by a 
factor of ~ 10 between apogalacticon and perigalacticon, there 
are short spikes coinciding with disk crossings, when the tidal 
field is compressive and > 10 times stronger than otherwise. 
The duration of these shocks is given by the disk scale height 
divided by the cluster's vertical velocity and amounts to lOMyr 
or less. The strongest shocks, which occur at disk crossings near 
perigalacticon, happen roughly once per orbit. Actually, the last 
such shock occurred 146 Myr ago, while the next one is due in 
110 Myr. 

With numbers collected in the previous subsections, we find 
from equation 10 that a typical relative change of kinetic en- 
ergy amounts up to ~ 30% for a disk shock near perigalacti- 
con, while equation gives i s h ~ 2 Gyr for the characteristic 
shock-evolution time scale. 

Another quantity of interest is the internal dynamical time, 
which we may estimate from the velocity dispersion to currently 
be ~ 80 Myr at the half-mass radius = 30 pc. After a strong 
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Fig. 3. — Strength of the tidal field along the orbit of Pal 5 (Fig. Q as mea- 
sured by the eigenvalues of c9 2 <J> / dxidxj . Positive (solid) and negative (dot- 
ted) eigenvalues correspond to compressive and stretching tidal forces, respec- 
tively. The little arrows on the top and bottom of the plot indicate moments 
of disk crossing and pericentric passages, respectively. 4 = corresponds to 
today. 

Table 1 
Time Scales for Pal 5 



dynamical process 




time scale 


duration of shocks 




< lOMyr 


internal dynamics 


^dyn 


~ 80Myr 


time between shocks 


idisk 


«3()0Myr 


shock-driven evolution 


*sh 


~ 2 Gyr 


two-body relaxation 


tih 


~ 20 Gyr 



NOTE. — Cluster internal time scales are esti- 
mated at the half-mass radius r\ ~ 30 pc. 



shocks and effects driven by two-body relaxation (partly caused 
by close encounters) are much less important (stellar mass loss 
affects clusters only in the first Gyr of their life). 

3.1. The Galactic Gravitational Field 
In our simulations, we considered only one choice for the Galac- 
tic gravitational potential; any effects due to variations of the 
Galactic potential are beyond the scope of this paper, but see 
the discussion in section l7"3l Because the evolution of the sim- 
ulated cluster is dominated by disk shocks, it is important to 
model the Galactic disk as realistical ly as possible. We used 
model 2 of iDehnen & Binnevl {l998), which contains a com- 
pound of three exponential disks for the thick and thin stellar 
disk and for the ISM, respectively, and two spheroidal compo- 
nents for the bulge and halo. The potential parameters have 
been determined to fit all observational constraints known in 
1998 and has a disk scale length of 2.4kpc, consistent with re- 
cent studies of the infrared background light distribution in the 
Milkv Wav (iDrimmel & Spergell200ll) . 

3.2. Initial Conditions 
To create sui t able initial conditions, we used King m odels 
( lMichielll963l iMichie & Bodenh eimerlll963t iKindl 19661) with 
baryce nter at the position and velocity of the orbit of Pal 5 (see 
section l2~Tl and Fig.Q ten radial orbital periods (corresponding 
to 2.95 Gyr) ago. King models have three free parameters, two 
scales (size and mass) and a shape parameter Wq, which is a di- 
mensionless measure for the depth of the gravitational potential 
and equivalent to the concentration c. As independent parame- 
ters, we used Wq, the total initial mass Mq, and the galactocen- 
tric radius 



disk shock, the cluster needs a few dynamical times to settle into 
a new equilibrium. Since the dynamical time is not substantially 
shorter than the time between subsequent shocks we anticipate 
that the cluster hardly ever is in dynamical equilibrium. 

Yet another important time scale is the two-body relaxation 
time t rn , whic h we can estimate u sing Spitzer & Hart's (1971) 
formula (Binn ev & Tremainel 19871 eq. 8.71). At the half-mass 
radius, we find 



Rt 



' lim 



GMx 







(6) 



a 


3 




-1 




0.22 kms" 1 




_O.OlM pc -3 _ 




_O.3M _ 



t lh ~ 20 Gyr 

0.22kms _i O.O1A40PC - " 0.3Mq 

(5) 

which is too long for relaxation related processes to currently 
play an important role and implies j3 ~ 10 for Gnedin et al.'s 
shock-importance parameter. 

We have summarized these various time scales in Table ^ 
Note that the estimates for internal time scales (td yn , ^sh, ^rh) 
must be regarded as crude ones, because they are based on the 
assumption of dynamical equilibrium, which we just saw is not 
really satisfied. We conclude, nonetheless, that disk shocks for 
Pal 5 are fast (shock duration < tdyn) and frequent (tdyn 
idisk) and entirely dominate the evolution of this object (t s h < 
tad. 

3. THE SIMULATIONS 

In order to model the tidal disruption of low-concentration glob- 
ular clusters in general and Pal 5 in particular, we performed 
collision-less iV-body simulations, i.e. with softened gravity, 
suppressing short-range interactions. This is entirely justified, 
because the evolution of these systems is dominated by disk 



with v c — 220 kms -1 , which means that the limiting radius 7 
Him of the model equals its tidal radius rtid (equation^ when 
at a distance R t from the Galactic center. With the choice of R t 
as parameter (rather than, say, ru m ) the importance of the tidal 
field is largely independent of the cluster mass when the other 
two parameters are kept fixed. 

We considered 1056 models from a grid of points in the 3D 
parameter space. For Rt, we took twelve values between 7.5 
and 13kpc at steps of 0.5 kpc. These are all larger than the 
perigalactic radius of Pal 5, guaranteeing that tidal interactions 
are important. For Wq, we used eleven values between 1.75 
and 4.25 in steps of 0.25, corresponding to low concentrations 
c between 0.46 and 0.88. Finally, for the initial mass the eight 
values 8, 10, 12, 16, 20, 24, 28, and 32 thousand Solar masses 
have been employed. 

Clearly, these initial conditions are, as often with iV-body 
simulations, somewhat artificial. We do not mean that 3 Gyr ago 
Pal 5 was well fit by a King model, rather we expect that it might 
already have undergone severe tidal disturbances. However, the 
tidal field is likely to quickly erase any specialities of our initial 
models, except, of course, the lack of possible tidal tails origi- 
nating from earlier epochs. Moreover, even if one were willing 
to improve on the initial conditions by integrating for the whole 
lifetime of a cluster, one had to know the detailed history of 



7 The outermost radius of a King model is often called its 'tidal radius' with 
the idea that a globular cluster ought to be tidally limited. However, as Pal 5 
clearly demonstrates, the limiting radi us m ay be different, actually larger, than 
the theoretical tidal radius, see section l231 



Modeling the disruption of the globular cluster Pal 5 Galactic tides 



5 



tj [km s l ~\ 



GM/ct 2 [pc] 



J? t =8 kpc 

C = 4 



fl t =9 kpc 
W„ = 3.5 



i? t =10 kpc 



i? t =ll kpc 
W=2.5 



R t =12 kpc 
W„=2 




-2.5 -2 -1.5 -1 -0.5 
t [Gyr] 



2 -1.5 -1 
t [Gyr] 



-0.5 



Fig. 4. — Time evolution of the mass (left), velocity dispersion (middle), and virial radius (right) of the simulated globular cluster as function of time for five represen- 
tative models with parameter Rt and Wo as indicated and Mo = 12000 Mq . The dotted vertical lines indicate disk crossings where the the tidal force is exceptionally 
large, quantified by the absolute largest eigenvalue of / dxidxj exceeding 7000 (km s — 1 kpc - 1 ) 2 , see Fig.|5] The two simulations with (R t , Wo) = (11, 2.5) 
and (12, 2) have been stopped before t = 0, because the number of bodies in the cluster dropped below 1000. 



the Galactic tidal field and would require much more CPU time 
without gaining much scientific significance. 

3.3. Technical Details 
We used the publicly available iV-body code gyrfalcON, which 
is based on Dehnen's (2000, 2002) force solver falcON, a tree 
code with mutual cell-cell interactions and complexity O(N). 
falcON not only conserves momentum exactly, but also is about 
10 times faster than an optimally coded Barnes & Hut (1986) 
tree code. 

Each simulation used TV=32000 bodies and a softening 
length of e = 3pc with what we call the P2 softening kernel, 
i.e. the Newtonian Greens function $ = — G/r was replaced by 



$(r) = - 



G 



s/r 2 + e 2 



1 



3e 4 



2 ( r 2 +e 2) ' 4( r 2 +e 2)2 



The density of this kernel falls off more steeply at large r than 
that of the standard Plummer kernel and becomes actually neg- 
ative (though with absolute value smaller than for the Plummer 
kernel) such that force bias is substantially reduced, see also 
Dehnen (2001). With these softening parameters, the maximum 
possible force between two bodies is equal to that for Plum- 
mer softening with e w 1.4pc, which hence might be called the 
'equivalent Plummer softening length'. 

The time integration was performed either for 2.95 Gyr, i.e. 
until today, or until cluster destruction. In practice, a cluster 
was considered destroyed if the number of bodies within one 
initial limiting radius from the cluster center dropped below 
1000 ~ 3% of the initial number. We used the leap-frog in- 
tegrator with step size of 2~ n Gyr rj 0.5 Myr and a block-step 



scheme that allowed up to 8 times smaller steps. The individual 
step sizes r where adjusted in an almost time-symmetric way 
such that on average r = a _1 kpcGyr _1 , where a denotes the 
modulus of the acceleration. This ensured that disk crossings, 
which have a duration of only a few Myr, are accurately inte- 
grated. With these settings, the energy of an isolated cluster 
was conserved to 0.8% over the period of 2.95 Gyr. One full 
simulation corresponds to about 6000 block steps and required 
about 75 min of CPU time on a linux PC (AMD, 1800MHz). 
1190 hours of CPU time were spent on the whole set of 1056 
simulations, 361 of which were halted before t = 0, because 
the cluster was found to be dissolved. 

3.4. Cluster Evolution in the Tidal Field 
For five representative simulations, ranging from medium con- 
centration and small extent to low concentration and large ex- 
tent, i.e. from least to most vulnerable to Galactic tides, Figure^] 
shows the time evolution of the cluster mass M, velocity disper- 
sion <7, and the ratio GM/a 2 . They have been computed from 
those stars that are within the original limiting radius from the 
cluster center, which was determined iteratively as barycenter 
of the same stars starting the iteration with the set of stars from 
the previous time step. 

3.4.1. The Mechanics of Tidal Disk Shocks 
As is evident from this Figure, the cluster evolution is driven 
by disk shocks (the instants of the strongest of which are in- 
dicated by dotted vertical lines), which may actually quickly 
destroy the cluster, depending on its initial state. Each strong 
disk shock causes an almost instantaneous increase of the clus- 
ter's velocity dispersion (middle panels of Fig. |4j by 8-25%. 
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Fig. 5. — Relative mass-loss and shrinking rates plotted versus each other (top 
left) and versus Wo for model clusters that have not been dissolved after 2 Gyr 
(for models dissolved earlier the uncertainties in the shrinking rate can be con- 
siderable). M and r v i r have been obtained by straight-line fits over the entire 
time interval, excepting the 0.5 Gyr before cluster dissolution, if applicable. For 
clarity, only models with Mo = 24000 Mq are shown (other values for Mo 
yield very similar results). The symbols refer to different values for Rt as indi- 
cated. 

This corresponds to an increase of the cluster kinetic energy by 
16-50% and pushes the cluster out of virial equilibrium. Stars 
which have been accelerated beyond their escape velocity are 
lost from the cluster, resulting in a drop of the cluster mass (left 
panels of Fig. [4}, which is delayed from the instant of the shock 
by the time needed for the stars to drift out of the cluster. 

In response to the heating by the disk shock, the cluster also 
expands. This together with the loss of the fastest stars re- 
sults in a substantial drop of the velocity dispersion by approx- 
imately twice the amount of the initial increase. This drop oc- 
curs roughly on a dynamical time scale, i.e. slower for the less 
concentrated and/or more extended clusters. Subsequently, the 
velocity dispersion shows damped oscillations until the cluster 
settles into a new equilibrium with velocity dispersion lower 
than before the shock. The damping is due to the fact that 
stars oscillate at different frequencies and is weakest for low- 
concentration clusters, since for those the range in orbital fre- 
quencies is smallest. For clusters of low concentration and/or 
large extent, the settling into a new equilibrium is forestalled 
by the next disk shock, i.e. these clusters are never in a state of 
dynamical equilibrium. 

3.4.2. Mass Loss and Time Evolution 
A remarkable observation from Fig.|4]is the fact that the orbit- 
averaged cluster mass decreases almost linearly with time over 
the durations simulated 8 . This implies that, at least for low- 
concentration clusters, disk-shocking induced mass loss quickly 
destroys the cluster. The orbit-averaged mass-loss rate depends 
strongly on the size, but also on the concentration of the clus- 
ter, spanning one order of magnitude in M/Mo amongst our 
simulations, see the bottom left panel of Fig. [5] 



8 We found the second time derivate M to be small in the sense that 
| M I M | > 7 Gyr with no preference for a sign . 



The linear mass loss is accompanied by an average decrease 
in velocity dispersion. In the right panels of Fig.|4] we plot the 
virial radius 

r vir = GM/a 2 (7) 

of the stars within the initial limiting radius of the cluster 9 . For 
most simulated clusters that survive until the present time, we 
find a modest decrease of the virial radius. In the top right panel 
of Figure[5] we plot the relative shrinking rate, — rVir/rvir(i = 
0), versus Wo for various values of R t but fixed Mq (the de- 
pendence on Mo is very weak at fixed R t and Wo). Most model 
clusters shrink slightly in response to the disk shocks with a rate 
that is primarily a function of Wo, or, equivalently, of the initial 
cluster concentration. Only very low-concentration clusters do 
not shrink but may even expand, see also the lower two simula- 
tions in Fig.0] The top left panel of Figure |5]plots the relative 
rates versus each other. While there is a clear tendency for mass 
loss to anti-correlate with shrinking, the spread in this diagram 
is considerable. 

We like to point out that the shrinking rates are very small, 
of the order of a few percent per Gyr, too small to protect the 
cluster from continued tidal stripping, i.e. self-limitation of disk 
shocks does not apply here, as is also evident from the un- 
damped mass-loss rates. 

3.5. Tidal Tail Morphology 
The stars lost from the cluster form two tidal tails, one leading 
the cluster and one trailing it. We will discuss the dynamics of 
the tidal tail in detail in the next section. Here, we want to give 
an impression of the morphology of these tails. 

For the model whose remnant cluster best fits Pal 5 (see sec- 
tion ^. 1> . hereafter 'model A', Figure|6]shows snapshots of (one 
out of ten) stars in tail and cluster: the particle distribution pro- 
jected onto the (instantaneous) orbital plane at time intervals of 
125 Myr. At the first time shown, the cluster has not yet been 
tidally shocked and all the particles are still inside the cluster, 
after that, the two tails develop. 

There are several notable observations possible from this 
Figure. First, near the cluster, the tail shows the typical S-shape, 
also visible in the observations of Pal 5, which originates from 
the trailing tail being at larger and the leading tail at smaller 
distance from the Galactic center. 

Second, the tail morphology varies between thin and long 
near perigalacticon (e.g. t = —0.45 Gyr), and short and thick 
near apogalacticon (e.g. t = —2.075, —0.075 Gyr), which is 
due to variations of the orbital velocity and hence the separation 
between tail stars and cluster. This is analogous to cars on a 
motorway: slow driving with short separations as in a traffic jam 
and fast driving with large distances correspond to the situations 
near apogalacticon and perigalacticon, respectively. 

Third, the tail neatly huddles against the cluster orbit at most 
of the times; only rarely near apogalacticon do some stars of 
the trailing (leading) tail appear closer (further). Fourth, occa- 
sionally near apogalacticon, the tidal tail shows a streaky struc- 
ture. Each of these streaks corresponds to a swarm of stars that 
has been set loose by the same disk shock and fills apparently 
a streak-like phase-space volume, which occasionally projects 
also onto a streak in configuration space shown in the Figure. 



9 One might instead want to use only stars that are bound to the cluster in 
the presence of the tidal field. However, due to strong changes of the tidal field 
during shocks, a star may switch from bound to unbound and back, resulting in 
spurious variations of the estimated quantities. Moreover, the contamination by 
unbound stars is small and does not significantly affect our estimates. 




Fig. 6. — Morphology of the simulated tidal tail in the simulation with parameters Rt = lOkpc, Wo = 2.75, and Mo = 20000 Mq (also called 'model A') at 
various times, given in Gyr in each panel. For every tenth body, we plot the position (in kpc) projected onto the instantaneous orbital plane of the cluster, whose 
motion is towards the right in each panel. The thin line indicates the cluster orbit. 



4. THE DYNAMICS OF TIDAL TAILS 

The gravitational influence of the cluster ceases not abruptly 
beyond its tidal radius, but it is completely negligible for stars 
that have drifted away from the cluster by more than a few tidal 
radii. Similarly, the self-gravity within the tidal tail is unimpor- 
tant, simply because the tidal-tail stars move too fast w.r.t. each 
other 10 . Thus, to very good approximation, the motion of a tidal 
tail star is governed by the gravitational potential <f> of the Milky 
Way only. 

For a time-invariant and smooth Galactic gravitational field, 
such as adopted in the simulations (but see section l73l for a dis- 
cussion of deviations from this ideal), a tidal-tail star moves on 
a Galactic orbit which is very similar to that of the cluster itself. 
The o ffset in orbital energy may be estimated to be ( Johnston 
1998) 

<5£ orb ~±r tid |V$|. (8) 

This follows from the fact that when the star leaves the cluster 
its cluster-internal energy approximately vanishes and its ve- 
locity is very similar to that of the cluster, which leaves only 
the difference in potential energy as significant contribution to 
SE oth . 

In reality, the star may leave the cluster with a higher than 
just escape energy, and |&E or b| follows a distribution with a tail 
to larger values, but in any case IJE'orbl <C |i?orb|- Thus, a star 
in the tidal tail moves on an orbit with slightly higher (lower) or- 
bital energy than the cluster, and hence slightly longer (shorter) 



1[, This may be different in tidal tails emerging from galaxy interactions, 
where self-gravity in the tails, possibly in conjunction with gas-dynamics, may 
lead to the for mation of b ound clumps that then become so-called 'tidal dwarf 
galaxies' (e.g. Due et al. 2000). 



orbital period. Over time, the longer (shorter) orbital period re- 
sults in a lag (lead) of the star w.r.t. the cluster, i.e. positive and 
negative SE OI b result in a trailing and leading tidal tail, respec- 
tively. Since the trailing tail moves on a higher orbital energy, it 
has a larger radius. 

If we assume that a tail star leaves the cluster at a larger or 
smaller galactocentric radius, but the same Galactic azimuth and 
latitude, and the same velocity as the cluster, then we expect it 
to have the same orbital eccentricity as the cluster and to move 
in the same orbital plane (paper II). In reality, we expect devia- 
tions from this ideal, but, of course, the orbital eccentricity and 
inclination should only slightly differ from that of the cluster 
orbit. 

4.1. Coordinates for Tidal-Tail Stars 
Since the orbits of a tidal-tail star and of the cluster are so sim- 
ilar, the position of a tail star may be referred to the cluster's 
trajectory. In paper II, we gave a simple model for the kine- 
matics of stars in the tail, which was based on the assumption 
that its orbit has the same eccentricity as that of the cluster and 
differs only in energy. If we further assume that the Milky Way 
has a flat rotation curve, we have 



\Xorb\ 



\SP\ 



(9) 



where SP is the difference between the orbital period of the tail 
star and the period P of the cluster, while Sx = x — cc or b is 
the positional offset of a tail star from the cluster trajectory at 
the same orbital phase, not the same time, see Fig. for an 
illustration. 

Since at the same phase Sx || x OI ^, equation l|9} gives a sim- 
ple relation between the perpendicular positional offset of the 
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Galactic center 

Fig. 7. — Illustration for the definition of tail coordinates, see section l4*Tl for 
details. 

tail star from the cluster trajectory and the period difference, 
which is directly related to the positional offset along the tra- 
jectory. In paper II, we used this simple relation to estimate the 
mean drift rate, and consequently the mass-loss rate, from the 
mean radial orbital offset of the tidal tail of Pal 5 from its orbit. 
For a flat-rotation-curve Galaxy, we further have (for SP <C P) 



SP 
P 



= cxp 



SE n 



- 1 



SE 



(10) 



Motivated by this model, we now introduce simple coordi- 
nates for stars in the tidal-tail. For a tail star at position x, we 
first find the reference point x or \> on the local cluster orbit by 
minimizing \x x a; or b(r)| w.r.t. the time parameter r along the 
cluster orbit, see Figure for an illustration. With r defined 
such that a; or b(0) = xq, we can use the time r directly as the 
tail star's (relative) coordinate along the tail. For the trailing 
(sketched in Fig. [7) and leading stars, r < and r > 0, re- 
spectively. Note that, in contrast to the true positional offset 
between cluster and tail stars, r is not subject to 'seasonal' vari- 
ations between apogalacticon and perigalacticon seen in Fig. [6] 
since these are caused by variations in orbital velocity. 

As further coordinates we use the components 



fe|| = Sx ■ aj<jrb/|a5orb|) 
Sx± = Sx ■ iorb/l-E'orbl 



(11a) 
(lib) 



and C_L ~ 0. That is, not only should these coordinates be 
conserved, i.e. all the kinematics is contained in their construc- 
tion, but they are related in a very simple way. However, the 
model on which equation (|9jl is based ignores differences in or- 
bital eccentricity and inclination between tail-star and cluster. 
Since the orbital period is mainly a function of orbital energy, 
whereas the orbital position at fixed phase also depends sensi- 
tively on the eccentricity and inclination, we expect deviations 
from our ideal to mainly affect dl3at and (± w 0, but not so 
much (tllbl . 

4.2. Kinematics of the Tidal Tail 
We are now going to investigate in some detail the distribu- 
tion of tidal-tail stars over the relative coordinates (77, (m, 
(5£orb, and t\, thereby also testing the simple relations Jl 3I > 
above. For this purpose, we concentrate on the one simulation 
already used in Fig. |6] above, which is one that gives a good 
description of Pal 5 (see Section [5j. In this model the cluster 
originally had ru m = 53 pc. 

In the top left panel of Figure [8] we plot the dimensionless 
drift rate 77, computed at the end of the simulation, vs. the time 
t\ of loss from the cluster (evaluated as the last instant when 
the star crossed r\i m outwards). Clearly, each disk shock trig- 
gers the loss of a swarm of stars with a spread in r\ which is 
almost symmetric for trailing and leading tail. The morphol- 
ogy in (77, t\) of each of these shocks is the same. First the stars 
with largest \t)\ escape, because they have highest velocity, and 
somewhat later those with smaller \t)\. The smallest |ry| of any 
tail star decreases towards later times, because the cluster loses 
mass and escaping from it becomes easier. 

The drift rates 77 follow a broad distribution with a factor of 
four between minimum and maximum. Therefore, a swarm of 
tail stars set loose by one disk shock quickly disperses along the 
tail. In particular, the fastest drifting tail stars will soon run into 
the slowest drifting ones lost at the previous shock. This implies 
that at any distance from the cluster the tail is a superposition of 
stars lost at various epochs, illustrated also in Fig.llOlbelow. We 
may estimate the number of shocks contributing to the tidal tail 
at some 'distance' r by dividing the difference between max- 
imum and minimum t — t\ by the average time idisk between 
shocks 



of Sx parallel to x or b (in the orbital plane) and parallel to the 
angular momentum vector L rb = ^orb X iorb (perpendicular 
to the orbital plane and the paper in Figs. [6] and 0. 

We may go a step further and define the dimensionless coor- 
dinates 



77 = r/it-k), 
C|| = Sx\i/\x olh \, 

C_L = Sx±/\x OT b\, 



(12a) 
(12b) 
(12c) 



where t\ is the time when the star was lost from the cluster. 
77 essentially is a dimensionless drift rate, while £11 and fj_ are 
dimensionless offsets in and out of the orbital plane. 

The time t corresponds to a phase difference S6 = tuj, 
which originates from the frequency difference Sco acting over 
the time t — t\, i.e. 59 = (t — t\)Sui. We thus have 77 = Suj/uj = 
—SP/ P and equations (|9jl and dlOt yield 



C11 



-77 

SE n '~ 2 



^orb 



K 



(13a) 
(13b) 



V 'mm 



^disk 4?7 

min L d.isk 



(14) 



With 77 m i n ~ 0.005 and idisk ~ 300 Myr, this gives n ~ 

0. 5(r/Myr). For the current situation of Pal 5, this translates 
into n ~ 2(Z/deg) with angular separation I from the cluster 
(assuming ?; or b ~ 100 km s -1 for the cluster's orbital velocity). 

The maximum drift rate is still rather low: about 2 percent, 

1. e. it would take ~ 50 orbits ~ 15 Gyr for the fastest first es- 
capees to cover a full orbital period, i.e. to fill a great circle on 
the sky. 

The top-right panel of Fig. [8] shows that there is indeed a 
tight relation between <5-E rb and 77, as expected from our equa- 
tion dl3b> . However, there is a considerable spread, in particular 
at small \r}\, which persists if we restrict the analysis to larger 
|£i| (excluding stars still in the neighbourhood of the cluster). 
There is also a systematic deviation in the sense that the lin- 
ear relations for the trailing and leading tail have non-vanishing 
zero points. 

The bottom panels of Fig. [8] show histograms of t\, i.e. the 
mass-loss rate, and <5i? rb- A comparison with Fig.[5]shows that 
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Fig. 8. — Distribution of tidal-tail stars over the dimensionless drift rate r\, time t\ of loss from the cluster, and orbital energy offset <5-E or b from the cluster. The values 
are obtained at t = from model A (Rt = 10 kpc, Wo = 2.75, and Mq = 20000 Mq). Tidal-tail stars with t\ > —0.2 Gyr have been omitted. The thin lower 
histogram in the bottom right panel is for stars with t\ £ [—2.5,-2.4] Gyr. In the upper right panel only one out of four bodies is plotted. 



the amount of mass lost in a shock correlates with the strength 
of the tidal field, as expected. The distribution in 6E OI b is, of 
course, a superposition of the distributions created by each of 
the disk shocks and the thin histogram shows the distribution 
from the single shock at —2.45 Gyr. The fact that the inner 
edge of the distribution (towards small |<5-E rb|) is not abrupt 
is mainly due to the fact that the cluster weakened with time, al- 
lowing smaller |<5.E or b| at later epochs. The typical scale for 
dEorb is ~ 300 (km s— 1)2. This value indeed follows from 
equation (jSJi for a cluster with rtid ~ 60 pc at R ~ 10 kpc in 
a logarithmic potential with v c ~ 220 km s -1 . The distribution 
in SE OI b is not quite, but still remarkably symmetric. 

In Figure|9j we plot the dimensionless orbital offsets (ji and 
(± vs. r\ computed at two different epochs, once near apogalac- 
ticon and once near perigalacticon. Quite obviously, the relation 
T) w — Q J 1 3at does not generally hold, only for the trailing tail 
at t — we find reasonable agreement. We also find that (± 
deviates from zero substantially in the sense that typical val- 
ues are significant compared to rj. A comparison between the 
two epochs clearly shows that 0i and (± are in fact not con- 
served. This, together with the violation of Q w — J7and£j_ sa 
implies, in accordance with the discussion following equations 
i!3l . that the orbital eccentricities and inclinations of tail stars 
differ significantly and systematically from that of the cluster. 
On the other hand, 6E or ^ shows only marginal changes with 
time (not shown in a Figure), while rj changes significantly with 
time, but still much less than 100%, in contrast to 0i and (±. 

We may thus conclude that in practice the dimensionless drift 
rate r\ is reasonably well conserved and correlates well with 
<5i? or b, while Qi and (± are subject to substantial time evolu- 
tion. These conclusions are supported by all simulations, not 
just the one presented in Figures|8]and|9] 

4.3. Stellar Density along the Tail 
The linear stellar density q — dM/ds along the tail, or, simi- 
larly, the circularly averaged surface density E(i2) of stars (tail 
& cluster) are easily observable from Pal 5 (see Fig. [2] paper II). 
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Fig. 9. — Plots of fx (top) and Qi (bottom) vs. r) for model A (the same as 
in Fig. [8} at t = {left, near apocenter) and t = — 0.45 Gyr (right, near 
apocenter). The thin lines indicate the relations 1 1 3al and (j_ = 0, which 
are expected from naive modeling. For clarity, only one out of four bodies is 
plotted. 



A study bv lJohnston. Sigurdsson & Hernauis^ (fr999) of the for- 
mation of tidal tails due to Galactic tides acting on globular clus- 
ters or satellite galaxies suggest that S oc R , or equivalently 
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Fig. 10. — Distribution of tidal-tail stars in t and time t\ of loss from the cluster 
for model A (the same as in Fig.|8}. The two vertical lines bracket the range in 
t inaccessible to stars with t\ < —2.95 Gyr. That is, stars lost prior to the start 
of the simulation, which are, of course, not contained in our model, contribute 
only to [t| > 19Myr. 

q w const 11 

We may identify three factors that determine g(s): the mass- 
loss rate M; the distribution f(rf) of dimensionless drift rates, 
which may be a function of the time t\ of loss from the cluster; 
and the orbital kinematics that map the phase-distance r into 
physical distance s along the orbit. The first two factors, which 
depend mainly on the general properties of the orbit (period and 
perigalactic radius) and the cluster, determine the distribution of 
tail stars over r 



dtjMfa) / 6rif(ri;k)S{T-ri(t-ti)). (15) 



AM 
17 



The last factor depends on the current position on the orbit and 
yields 

, , 1 dM 

Q(s;t) = - --r— (16) 

Krb(i")| dr 

where s — f Q |u or b ("?"') I dr' calculated at the current time t. 

A constant g(s), such as found in the simulations by John- 
ston et al., emerges from a constant (orbit-averaged) M in con- 
junction with a time-invariant distribution /(ry) of drift rates and 
\votb\ ~ const over the extent of the tail. This latter condition 
is not satisfied whenever the tidal tail spans a significant range 
in phases on an eccentric orbit. 

As we saw in Fig. [8] the distribution over drift rates r\ is ac- 
tually not time invariant for our models: the drift rates decrease, 
because the cluster is significantly weakened by the continued 
loss of stars. With a constant (orbit averaged) mass-loss rate, 
this results in dM/dr slightly decreasing with |t|. This is il- 
lustrated in Figure which plots the distribution of tail stars 
over t and t\. In the top panel, it is evident that the typical 
|t| increases slightly faster than linear with t — t\ (the inner en- 
velopes are not straight lines). The plot of dM/ dr in the bottom 
panel shows (i) some 'bumpiness' at the < 20% level and (ii) a 
slight decrease with |t| (data for |t| > 19Myr are incomplete 



11 In another study Johnston, Choi & Guhathakurta 1 2002 ) reported much 
steeper profiles (S oc R~ a ). However, these occurred in much more massive 
satellites and over larger parts of the orbit and are not relevant to the situation 
studied here. 
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Fig. 11. — For all 695 simulations with a cluster surviving until t = 0, we 
plot against each other: the average offset in orbital energy <5-E or b between 
tail stars and cluster; the average absolute dimensionless drift rate r\\ and the 
weighted mean of the initial and final mass of the cluster. SE OI ^, and r\ have 
been calculated at t = 0. The thin curves in the two left panels follow a power- 
law with exponent 1/3. 

since tail stars with t\ < —2.95 Gyr are lacking from our mod- 
els, see also the discussion at the end of section EOl i. For the 
linear density g, the latter implies a slight decrease with s and a 
steeper than R~ x decrease for E. Thus, we expect non-constant 
g(s) even for a (orbit-averaged) constant mass-loss rate and a 
tail spanning a small range of |u rb|. if the mass loss was sub- 
stantial. 

For clusters or satellites which undergo only non-substantial 
mass loss, the distribution over r] should indeed be independent 
of time and consequently dM/dr w const, consistent with the 
results of the aforementioned study by Johnston et al.. 

The bumpiness of dM/dr, which of course translates in a 
bumpiness of g(s), is caused by individual shocks and the fact 
that for small r only few of them contribute. These bumps are 
expected to be symmetric, i.e. occur in trailing and leading tail 
at about the same distance from the cluster. For larger t or 
distances s from the cluster these inhomogeneities are averaged 
out by the superposition of many swarms of stars lost from the 
cluster at different shock events. 

4.4. Drift Rates and Tail Kinematics of the Models 
In order to investigate the dependence of the tail properties on 
the model parameters, we computed for each model of a cluster 
that survived until today the mean and dispersion of r\, Qi, 
and <5£orb evaluated at t = for both the leading and trailing 
tail. 

The estimate (|8j for i5i? rb suggests that (|<5-E rb|) depends 
only on the mass M of the cluster at the instant of the star 
leaving it. Because our models are subject to substantial mass 
loss, their mass is not conserved. However, as the bottom- 
left panel of Figure ITTI demonstrates, the average (|5S or b|) 
over all tail stars (lost over the whole integration interval) is 
strongly correlated with some time-averaged cluster mass. In 
fact (|<5-E rb|) oc M 1 / 3 , exactly as equations (|8} and Q predict. 
The constant of proportionality in this relation depends on the 
cluster orbit and the Galactic potential only, but not on the clus- 
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Fig. 12. — Distribution of the 695 simulated clusters over the fractions of stars 
within projected 16' that are also within projected 4' (x-axis) and 12' (y-axis). 
Note that 16' is the limiting radius of Pal 5 itself, whose data are shown as solid 
square. The open square marks the model which best fits the data for Pal 5 and 
is denoted 'model A' in the text. 

ter's intrinsic properties. The remaining two panels of Fig. II II 
show that the mean dimensional drift rate (|ry|) is proportional 
to M 1 / 3 , too, and to (|<5-E rb|)- The constant of proportionality 
in this latter relation is (196km/s)~ 2 , corresponding approxi- 
mately to v~ 2 , in accordance with equation dl3bl i. 

With this mass dependence, we obtain for the average density 
of tail stars over r from equation dl5> 



/dM\ M 



(17) 



with M denoting some measure of the cluster mass (final or 
time averaged). The constant of proportionality depends on the 
relative mass-loss rate M/Mq and hence is a function of R t and 
Wo, but hardly of Mq. Thus, at fixed (R t , Wq), a smaller cluster 
mass results in a larger fraction of tail stars at fixed r, i.e. phase 
offset from the cluster. 

5 . COMPARISON WITH PAL 5 

In this section we compare the 695 models with a surviving 
cluster to the observed properties of Pal 5 and its tail. The ob- 
jective is to investigate whether these data are at all consistent 
with our models and if so, with which type of models. 

5.1. The Cluster 
The most astonishing property of Pal 5, apart from the tails, is its 
size, which is about four times larger than expected for a cluster 
of its mass and on its orbit ( section l2~3l . Instead of comparing 
695 different projected surface density profiles, we plot in Fig- 
ure^]the fractions of the projected cluster mass within 16' (the 
limiting radius of Pal 5), that are also within the projected radii 
of 4' and 12'. Most of our models are in the upper-right corner 
of the figure: they have nearly 100% for both these fractions, 
meaning that most of their mass is contained within 4', the tidal 
radius of Pal 5 at perigalacticon. Thus, these simulated clusters 
have survived until today, because they are, in contrast to Pal 5, 
smaller than their tidal radii. 

There is, however, also a long tail of models with lower 
mass fractions within 4' and 12', corresponding to larger ex- 
tents. Some of these actually come close to Pal 5 in the sense 




Fig. 13. — Time evolution of mass (top) and velocity dispersion (bottom) of 
model A (with R t = 10 kpc, W = 2.75, and M = 20000 Mq). The dotted 
vertical lines are equivalent to those in Fig.|4] The cluster is destroyed by the 
first future disk crossing at t = 110 Myr. 
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Fig. 14. — Distribution of the 695 simulated clusters over the final projected 
mass within 16' and the line-of-sight velocity dispersion within 6'. Larger sym- 
bol sizes correspond to m odels that are closer to Pal 5 in their enclosed projected 
mass fractions tFig. 1121 : the largest symbol refers to model A. Note that for 
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that they have ~ 95% of their mass within 12' and ~ 60% 
within 4'. These models have initial parameters Wq ~ 2 — 3, 
Rt ~ (4.5 + 2W )kpc, and M ~ 20000 -Mq, while model 
clusters with R t larger than that (at fixed Wo) get destroyed and 
those with lower mass are even more diffuse or have been dis- 
rupted as well. Initially more concentrated models with Wo > 3 
and also R t ~ (4.5 + 2Wo) kpc form the tail in Fig.ll2ltowards 
smaller mass fractions; they have massive tidal tails but their 
clusters are not as extended as Pal 5. 

The model that best matches Pal 5 in Fig. ^]is that with 
initial parameters R t = 10 kpc, Wo = 2.75, and Mq = 
20000 M Q ; it is marked with an open square and will be de- 
noted 'model A'. This model already featured in Figures [6j [8] 
[9] and ^| and will be used in the remainder of the paper for 
illustrative purposes. 
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Fig. 15. — Top: distribution over distance r from the Sun and line-of-sight 
velocity for stars within projected 6' of the cluster center for the final state of 
model A (initial parameters: W = 2.75, R t = 10, and M = 20000M Q ). 
Bottom: distribution of the same stars over vi os and a Gaussian with cn os = 
0.38 km s" 1 (dotted). 

Figure [O] shows the time evolution of the cluster mass and 
velocity dispersion a of model A. All models which come close 
to Pal 5 in their projected mass fractions in Fig. 1121 show sub- 
stantial mass loss, which results in a drastic increase in the clus- 
ter's dynamical time idyn, reflected by the steady decline of 
a. The increase of t& yn slows down the cluster's response to 
disk shocks: the period of the (damped) oscillations that follow 
each shock (in Fig. [U] visible as oscillations of a subsequent 
to each dotted vertical line, see also the discussion in section 
13.4. \\ increases. Eventually, this period approximately equals 
twice the time between two shocks, namely those at t = —146 
and HOMyr, such that just in the middle between these two 
shocks, i.e. today, the cluster is in state of maximal expansion 
and minimal velocity dispersion. This property is common to 
all models with projected density profile similar to that of Pal 5. 

Thus, all those models that reproduce the overly large extent 
of Pal 5 also naturally, and causally connected, predict a very 
small velocity dispersion consistent with that actually observed 
for Pal 5. 

In Figure ^] we plot for all 695 surviving cluster models 
the final mass within projected 16' (the limiting radius of Pal 5) 
and velocity dispersion of stars within projected 6'. The mod- 
els with projected mass profiles most similar to Pal 5 (shown as 
larger symbols) form the lower envelope, i.e. have lowest line- 
of-sight velocity dispersion eri os at given projected mass. A di- 
rect comparison of these values for cri os with that determined for 
Pal 5 in paper I (see section l2~2l is problematic, as the latter has 
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Fig. 16. — Distribution of the 695 simulated tails with surviving clusters over 
(:r-axis) the ratio between mass in the tail out to projected 200' and the mass in 
the cluster (within 16') and (y-axis) the fraction of tail stars within 100' for the 
trailing (top) and leading (bottom) tail, respectively. The symbol sizes indicate 
how wel l the model matches the cumulative projected masses of the cluster 
(Fig. 1121 . The solid square represents the star-count data for Pal 5, which have 
a statistical uncertainty of about 0.05 in the y direction. The trailing tail of Pal 5 
has a substantial clump at 3 deg from the cluster, compromising the measures 
plotted, in particular that on the y axis. 

been obtained under the assumption that the line-of-sight ve- 
locity distribution was Gaussian (high-velocity tails have been 
modeled as binary contamination). 

In Figure we show for the model that best fits the pro- 
jected mass profile (the biggest symbol in Fig. 1141 the distribu- 
tion of stars within projected 6' of the cluster center over dis- 
tance r and line-of-sight velocity v\ os . The dotted line in the 
lower panel represents a Gaussian with eri os = 0.38 kms -1 . 
Clearly, the distribution has non-Gaussian wings, which orig- 
inate from stars that are just being lost from the cluster. Thus, 
the extra-Gaussian tails of the observed line-of-sight velocity 
distribution can at least partly be explained by cluster dynamics 
and not entirely by binaries, as we assumed in paper I. 

The very model shown in Fig.[2]not only gives the best fit 
to the shape of the Pal5's mass-profile, but also its total mass 
within 16' of 5177 Mq agrees nicely with the best estimate for 
Pal 5 of 5200 ± 700 M Q (sectionlO. 

5.2. The Tidal Tail 
5.2.1. The Radial Profile 
In Figure^] we compare the gross radial distribution of stars 
in the tail of Pal 5 with those of the 695 models with surviving 
clusters. On the cc-axis, we plot, for the trailing (top) and leading 
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Fig. 17. — Distribution of the 695 simulated tails with surviving clusters over (a>axis) the mean projected drift rate along the orbits and (y-axes) the mean (bottom) 
and dispersion (top) of the rectangular displacement p of the tail stars from the orbit. Only tail stars with distance (projected arc-length s along the tail) from the 
cluster within [24', 200'] have been used. The symbol sizes indicate how well the model matches the cumulative projected masses of the cluster (Fig. 1121 . See text 
for a comparison with Pal 5 



(bottom) tidal tail, the ratio of tail to cluster mass. Here, we 
rather arbitrarily truncated the tail at 200' to ensure that no stars 
are missing from the simulations due to our limited integration 
time (cf. the discussion following Fig. II 0i . On the y-axis of 
Fig- El we pl°t the fraction of the tail mass within 100'. This 
fraction is related to the density run along the tail: for a constant 
linear density, equivalent to S a r 1 with 7 = 1, we expect 
~ 0.46, while steeper density profiles with 7 > 1 yield larger 
values. 

We first notice a weak correlation between this ratio and the 
relative mass in the tail in the sense that (relatively) more mas- 
sive tails have higher fractions of stars inside 100', correspond- 
ing to larger 7. This is, of course, exactly what we expect from 
our discussion in £14. 31 

Regarding Pal 5, we see that while its leading tail falls into 
the locus of models, the trailing tail does not. This can be en- 
tirely attributed to the clump at ~ 3° from the cluster (paper II, 
see also Fig. [2). We also notice that those models which best 
fit the properties of the cluster itself (indicated by larger symbol 
sizes in Fig. I16i. have tails that are more massive by a factor 1.5 
- 2 and have a steeper radial density profile (but note that the 
statistical uncertainty on this latter property of Pal 5 is consid- 
erable). 

We should note here that for Pal 5, we actually cannot mea- 
sure mass ratios, but merely number count ratios, which are 
identical to the former only, if the stellar luminosity functions 
(LFs) in t he cluster and the v arious parts of the tail are the same. 
Actually, Koc h et al.l (12004) reported a significant lack of low- 
mass stars in the cluster as compared to its tails, which may 
explain (at least partly) the apparent discrepancy in the tail-to- 
cluster mass ratio. 

Thus, we conclude that simultaneous fitting of the projected 
radial profile of Pal 5 and its tails appears to be difficult, while 
each property on its own is inside the locus of the models. How- 



ever, statistical uncertainties and our lack of knowledge of mass 
(rather than number-count) ratios, i.e. of the LF, hampers more 
quantitative statements to be made. 

Another important point is that the dumpiness of the tidal 
tail of Pal 5, in particular the over-density at 3° in the trailing 
tail, is not reproduced in any of our models. We will discuss in 
section fOl the possible implication of this mismatch. 

5.2.2. The Offset from the Orbit 

In paper II, we measured the projected rectangular offset of tail 
stars from the projected (assumed) orbit. In order to compare 
our simulations to these measurements, we computed the same 
quantities for our models (after projection onto the sky as seen 
from the Sun) as follows. First, we determined the point on the 
projected orbit that minimizes the projected distance to the tail 
star. We then define p to be the projected distance with p > 
(p < 0) for tail stars at larger (smaller) Galactic latitude than 
the nearest point on the orbit. We define s to be the arc-length 
along the projected orbit with s > and s < for the leading 
and trailing tail, respectively. 

In Figure we plot, separately for the trailing (left) and 
leading (right) tails of the all the models with surviving cluster, 
the mean (bottom) and dispersion (top) of p versus the mean 
projected drift rate s/(t — t{). 

From this Figure, we see that, while the drift rates for the 
leading and trailing tail are of comparable size, the mean rect- 
angular distance \p\ is larger by almost a factor of two for the 
former. This is also evident from the projection of model A in 
the bottom panel of Figure The reason is partial intrinsic, 
i.e. due to the offsets £ji of the tail stars from the cluster orbit 
within the orbital plane (Fig. |9] left bottom panel), and partly 
due to projection effects (because the leading tail is closer to us 
and seen under a slightly smaller inclination angle). 

In contrast to the modeled tails, the values inferred in pa- 
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Fig. 18. — Projection of model A at t = (today) onto the sky in standard Galactic coordinates (bottom) and plots of the offset in distance modulus (middle) and 
radial velocity (top) of the stars from the cluster versus Galactic longitude. The thin line represents the orbit. Only every second body is plotted. 



per II for the mean rectangular offset of Pal5's tail stars are 
rather similar for the trailing and leading tail (11/.8 ± 0'.5 and 
— 10'. 1 ± 0'.8, respectively). However, this is essentially a con- 
sequence of the implicit assumption, made when adopting the 
orbit in paper II, that the rectangular offsets are similar in size 
between the two tails. The fact that the simulations consistently 
predict larger (\p\) for the leading than the trailing tail implies 
that the orbit of Pal 5 should actually be slightly different from 
the one used in paper II and here; compare, for instance, the 
bottom panel of Fig.^]with Fig. 9 of paper II. (Note, however, 
that the uncertainty in the distance to Pal 5 still dominates the 
uncertainty in the orbit). 

However, we may still compare the tails' perpendicular 
structure to the models by means of the dispersion in p. In 
paper II, we fitted a Gaussian to the distributions in p for the 
trailing and leading tails of Pal 5. When assuming that the dis- 
tribution is actually Gaussian, we obtain for the dispersions 
7'.8 ± 0'.5 and 7'.3 ± 0'.8 for the trailing and leading tail, respec- 
tively. Comparing these values with those in Fig.^J it appears 
that the tail of Pal 5 is at the upper end of the distribution. In 
particular, these values are somewhat larger than those for the 
models that best fit the cluster's cumulative mass distribution 
(big symbols in Fig. ll7t . 

The mean drift rate estimated for the trailing and leading tail 
in paper II was ~ 180 arcmin/Gyr, which is substantially more 
than predicted by the models, even by those with as high a dis- 
persion in p as observed for the tails of Pal 5. That estimate, 
however, was based on rather simplifying assumptions on the 
relation with p and, given the model results in Fig.^] a lower 
value for the mean drift rate seems more appropriate. In face 
of the uncertainties, we cannot make precise statements, except 
for saying that most likely the mean drift rate is in the range 
100 — 140 arcmin/Gyr, implying that the average tail star at the 
end of the observed trailing and leading tail was lost from the 
cluster about 3 ± 0.5 and 2.4 ± 0.4 Gyr ago, respectively. 



5.2.3. Clues for Further Observations of Pal 5 

In Figure [^J we plot the projection of model A onto the sky 
as well as the distribution of stars in their radial-velocity and 
distance-modulus offset from the cluster as function of offset in 
Galactic longitude. A comparison of the projection of model 
A with Fig. 9 of paper II shows, as already discussed above, 
that the orbit chosen in paper II is presumably not quite correct. 
This is important for any application where the orbit geometry 
is derived from the tidal tail. 

The other two panels of Figure^]plot the observable radial 
velocity and the distance modulus for the tail stars versus their 
galactic longitude. A measurement of either of these quantities 
for a sufficient number of tidal tail stars would constrain the 
acceleration in the Galactic halo at Pal5's present position and 
hence directly give the mass of the Milky Way within w 18 kpc, 
a quantity not very well known. 

From this Figure, we see that the distance-modulus offset of 
the tails is larger than that of the orbit, simply because of the 
generic shape of the tidal tails (see Fig. [6}. This enhances the 
signal when trying to measure this effect: over the currently ob- 
served arc-length of the tail (from about —6 to 4 in A(^cos6), 
see Fig. 9 of paper II) the distance modulus changes by about 
0.25 mag. The levelling off of the Am at I ~ 8° is because this 
is close to the apogalacticon of the orbit. This implies that in or- 
der to constrain the orbit of Pal 5, and hence the directly mass of 
the Milky Way, distance-modulus determinations should focus 
on^<8°. 

A similar effect is seen for the radial-velocity offset, which 
also shows a larger gradient than the orbit itself. The difference 
here is that the spread in Av r at any fixed position is larger, 
but also that gradient is maximal at apogalacticon, in contrast to 
Am whose gradient vanishes there. 

Both the gradient in Am and that in Av r are roughly propor- 
tional to the circular speed v c of the Galactic halo at the position 
of Pal 5. Thus, in order to determine v c to 10% accuracy, one 
would need to determine the (v r ) at ~ 5° on either side of the 
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cluster to about 4kms~ 1 accuracy. Note that more accurate 
individual measurements are not sensible, because of (1) the in- 
trinsic width of the tail and (2) possible binarity of the tail stars. 
For a determination of v c to 10% accuracy from the distance 
moduli at, say ~ 5° on either side of the cluster, an accuracy 
of 0.02 mag is sufficient. Both these accuracy requirements are 
not beyond current capabilities and a determination of the mass 
of our host galaxy within 18.5 kpc to 10 — 20% accuracy seems 
possible. 

6. DISCUSSION I: 
CLUSTER EVOLUTION DRIVEN BY DISK SHOCKS 

Our simulations were originally intended to understand the evo- 
lution and dynamical status of the globular cluster Palomar 5, 
but are also important for understanding disk-shocking driven 
globular-cluster evolution in general. 

The orbit of Pal 5 is rather eccentric with an apogalac- 
tic radius of about 19 kpc, near the cluster's present loca- 
tion. It carries the cluster to a perigalactic radius as low as 
5.5 kpc, which implies strong tidal shocks at disk-crossings 
near perigalacticon. For a low-concentration low-mass cluster 
like Pal 5, these shocks rather than internal processes (driven 
by two-body relaxation) dom inate the dynamical evolution 
ddnedin. Lee & OstrikeJl999h . 

A strong disk shock puts the cluster out of dynamical equi- 
librium by unbinding some stars and changing the (cluster- 
internal) energies of all the others. In addition to this energy 
mixing, there is also a general heating of the cluster. In low- 
concentration low-mass clusters, like Pal 5, this all happens 
much faster than their dynamical time, i.e. the shocks are im- 
pulsive. In response to this instantaneous heating, the cluster 
expands on its dynamical time scale, then contracts again, ex- 
pands and so on. These oscillations are damped, since the in- 
ternal dynamical times of the stars differ so that they get out of 
phase, and after a few dynamical times the cluster is in a new 
dynamical equilibrium. 

For very low-concentration clusters the internal dynamical 
times (1) are long and (2) do not differ much, so that the damp- 
ing is less efficient. This implies that for these systems the 
settling into a new equilibrium may be forestalled by the next 
strong disk shock. Hence, these clusters will never be in an 
equilibrium state, but are continuously rattled by the repeated 
tidal punches. 

6.1. Evolution of Extended Clusters 
We find that these shocks can be very efficient in driving the 
cluster evolution and may destroy the cluster within a few or- 
bital periods. The disruption is faster the lower the cluster's 
concentration and the larger its size compared to its perigalactic 
tidal radius. Thus, at smaller Galactic radii, the lifetime of low- 
concentration clusters is lower because of both stronger shocks 
and shorter orbital periods. 

We should emphasize here, that our models were not ini- 
tially, and never became, limited by their tidal radii. This im- 
plies that globular clusters extending beyond their theoretical 
(perigalactic) tidal radius and orbiting on eccentric orbits on 
which they experience strong disk shocks, will not quickly be 
'tidally stripped' down to their tidal limit. This result seems to 
be at odds with the general wisdom of globular cluster dynam- 
ics in general and with the results of Gnedin et al. in particular. 
These authors found, using Fokker-Planck simulations, that in 
response to disk shocking the clusters become more compact, 
which renders the shock-induced evolution self-limiting. 



However, we nonetheless think that our results are correct 
and not contradicting previous knowledge. Our simulations dif- 
fer in various ways from those of Gnedin et al., who consider 
disk shocks which are non-impulsive and rare (internal dynam- 
ical time much shorter than time between shocks). We attribute 
the difference in our findings mainly to the eccentricity of the 
orbit of Pal 5 as compared to that of NGC 6254 (the cluster stud- 
ied by Gnedin et al.). After being shock-heated near perigalac- 
ticon, clusters on eccentric orbits approach their new equilib- 
rium whilst being in the remote parts of their orbit, where the 
Galactic tidal field is much smaller (about ten times, see Fig. [3j 
than near perigalacticon, even away from the disk. As a con- 
sequence, the new equilibrium, which does not know anything 
about the conditions near perigalacticon, is not limited by the 
perigalactic tidal radius. 

This mechanism should not be very sensitive to the typical 
dynamical time of the cluster, since almost by definition of the 
tidal radius, the dynamical time of the outer parts of the cluster 
never is much smaller than the orbital time. 

We find that in response to the disk shocks the clusters shrink 
somewhat, except for very low-concentration clusters, which 
even expand, making them yet more vulnerable and acceler- 
ating their destruction. However, this shrinking is very mod- 
est, of the order of at most a few percent per orbit even for 
medium-concentration clusters, much too small to protect them 
from continued disk shocking. For clusters of higher concentra- 
tion and/or on less (or more) eccentric orbits, further studies are 
needed to investigate under which conditions tidal-force driven 
cluster evolution is self-limiting. 

One implication of this possibility of super-tidally limited 
globular clusters is that inferring the mass of the Milky Way 

from the limiting radii of globular cluster assuming they are 

l n l 1 1 ' ■ ■ ^ j 

tidally limited ( King 1962; Bellazzini 2004) is at best dangerous 
and may under-estimate the mass of the Milky Way. 

6.2. Implications for the Globular-Cluster System 

Many studies of the evolution of globular clusters and globular- 

cluster systems (e.g., [y^spe rini & Heggidl997H Baumgardt & Makinol 
2003; Fall & Zhang 2001) assumed that the sizes of globular 
clusters are initially upwardly limited by their (perigalactic) 
tidal radius. This assumption is based on the idea that the galac- 
tic tidal force field quickly results in such a tidal limitation, even 
if the clusters were originally larger. In view of our results, this 
assumption appears not generally justified, at least not for clus- 
ters on eccentric orbits. This implies that the above studies have 
probably underestimated the importance of disk shocks for the 
evolution of globular-cluster systems. 

Our simulations predict unambiguously that Pal 5 will not 
survive its next disk crossing in ~ HOMyr. If me make the 
assumption that there is a plausible initial condition that could 
have formed lOGyrs ago and is being disrupted on this orbit 
now, then the fact that we observe Pal 5 within the last per cent 
of its lifetime, suggests that we see only the last tip of a melting 
iceberg. The inner Milky Way may well have been populated 
with numerous low-concentration globular clusters, which by 
now have all been destroyed, except for those, like Pal 5, whose 
orbits have longer periods and carry them into the more remote 
parts of the Galaxy. Currently, we have little possibility of quan- 
titatively assessing the initial amount of low-concentration and 
low-mass globular clusters of the Milky Way galaxy. If tidal 
tails can survive for a Hubble time (but see section 1731 . we 
may be able to identify the debris of destroyed clusters as nar- 
row streams of stars, which in turn may be identifiable in deep 
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photometric, astrometric and/or spectroscopic surveys such as 
SDSS, RAVE and GAIA. 

6.3. Tidal Tails 
Each disk shock triggers the loss of a number of stars: those 
which have been accelerated beyond their escape velocity. 
These stars emerge with orbits whose energy is either slightly 
above or below that of the cluster. This offset in orbital energy 
is small compared to the orbital energy of the cluster itself (of 
the order of one per cent), so that the stars are not dispersed 
all over the Galaxy, but stay in a trailing and leading tidal tail, 
depending on the sign of the energy offset. This offset in orbital 
energy translates into an offset in orbital period and hence in a 
drift away from the cluster with a rate of the order of one per 
cent, i.e. a full wrap-around requires 0(100) periods. 

The orbital-energy offsets follow a broad distribution result- 
ing in a distribution of drift rates. This implies that stars which 
have been lost from the cluster at the same time will not stay 
together but disperse along the tail. At any distance from the 
cluster, the tail contains stars lost at various shock events. 

A constant drift rate in conjunction with a constant mass-loss 
rate results in a constant linear density of stars in the tidal tail 
or, equivalently, to a surface density £ oc r -7 with 7 = 1. This 
seems to be the generic property of tid al arms that results from 
contin ued tidal shocking iTohnston. Sigurdss on & Hern auist 
1999). However, when, as a result of continued mass loss, the 
cluster mass and hence the escape velocity drops considerably, 
stars with ever lower orbital energy offsets and hence lower drift 
rate can escape. In this case, the surface density exponent may 
be 7 > 1, in agreement with our empirical findings. 

7. DISCUSSION II: 
THE STATE AND FATE OF PAL 5 

7.1. Is Pal 5 in Equilibrium? 
Our simulations were primarily aimed at understanding Pal 5, 
predominantly because it is observed to possess a massive tidal 
tail. Apart from these tails, the most unusual property of Pal 5, 
which however has not been hitherto noticed, is its enormous 
size. The cluster is two times larger than its present theoretical 
tidal radius and even four times larger than its perigalactic tidal 
radius. This property can in fact be reproduced by (some of) our 
simulations, and can be understood in terms of the expansion 
that immediately follows each strong disk shock. All our model 
clusters that can reproduce the present size (and mass) of Pal 5 
are in a state of near-maximal expansion after the last disk shock 
which occurred ~ 146 Myr ago. 

Maximal expansion implies minimal kinetic energy, corre- 
sponding to a small velocity dispersion. Indeed, models that can 
reproduce cluster size and mass comparable to Pal 5 also have 
the smallest (line-of-sight) velocity dispersions of all models, 
namely a\ os ~ 0.45 km s -1 . We cannot produce models with 
a\ os much smaller than that, apparently because the tidal field 
limits the lowest velocity dispersions possible for a surviving 
cluster on any given orbit. 

The line-of-sight velocity distribution of the models is actu- 
ally not Gaussian, but contains some high-velocity tails result- 
ing from stars just leaving the cluster. Correcting for these, we 
find a Gaussian erio S ~ 0.38kms _1 for model A that best fits 
size and mass. This value is consistent with our observation for 
Pal 5 of cri os = 0.22+°;^ km s" 1 (paper I). 

All models that are comparable to Pal 5 in their cumulative 
mass distribution also predict that the cluster will not survive 



the next disk crossing in ~ 110 Myr. The models that best fit 
the cluster radial profile have a constant orbit-averaged mass- 
loss rate of ~ 5000 M Q Gyr _1 . Extrapolating this back in time 
implies an original mass of ~ 70000 Mq, though this may be 
an over-estimate, since the strength of the Galactic tidal field 
was probably weaker in the past. 

7.2. The Tidal Tail 
Our models also have massive tidal tails. However, models that 
best fit the properties of the cluster tend to have a fraction of 
stars in the tail within 200' that is about 50 — 100% higher than 
observed for Pal 5. Note however that in this comparison we es- 
sentially assumed implicitly that the cluster's luminosity func- 
tion equals tha t of the tail, which actually appears to be not quite 
correct llKoch et all20 04i. 

In our best fitting cluster models, ^55% of the tail stars 
within 200' are also within 100'. This corresponds to a surface 
density exponent of ~ 1.35 and is consistent with our observa- 
tions of Pal 5 (paper II). 

All our models have a rather smooth density distribution 
along the tail with structure of less than ~ 20% in amplitude. 
This is because, as mentioned in section l6~3l above. even though 
the stars are lost at discrete events, they quickly disperse along 
the tail, washing out the discreteness of their liberation from the 
cluster. 

This property of the model tails is in clear contrast to the 
observations of Pal 5, in particular its northern tail, which pos- 
sesses a substantial clump at about 3 deg and seems to fade away 
at ~ 6 deg, corresponding to a drift time of only ~ 3 Gyr (sec- 
tion !5.2.2> . 

7.3. What Causes Structure of the Tidal Tail? 
Why does the tidal tail of Pal 5 show this significant structure? 
Obviously, it must be caused by something not accounted for in 
our models. Evidently our models lack any Galactic substruc- 
ture, such as giant molecular clouds and spiral arms. Clearly, 
since the Galactic disk itself is rotating, different parts of the 
tidal tail encounter different substructure when crossing the 
disk. If the tail passes, say, a giant molecular cloud at a distance 
smaller than its length, different parts feel substantially different 
perturbations and the tail will develop some substructure. 

Another potential source of Galactic substructure that may 
affect the tidal tail are clumps in the dark-matter halo, which 
are actually predicted by stand ard ACDM cosmology (e.g. 
iMooreet alJll999t iKlvpin etal] ll999). or any super massive 
compa ct halo objects, such as proposed by lLacev & Ost riker 
( 1985). In fact, the very integrity of the tail of Pal 5 may impose 
useful upper limits on the amount, sizes, and masses of such 
possible dark substructure in the Milky Way halo. 

8. SUMMARY AND CONCLUSIONS 

In order to study the disk-shock driven evolution of low- 
concentration globular clusters in general and of Pal 5 in partic- 
ular, we performed 1056 iV-body simulations. All these have 
the same Galactic potential and orbit (that derived for Pal 5), but 
differ in the internal properties of the cluster: its initial concen- 
tration (measured by Wq for our King models), mass, and size 
(corresponding to our R t parameter). 

Our simulations demonstrate that disk-shocks near peri- 
galacticon are very efficient at driving the evolution of low- 
concentration globular clusters. In particular, the disk-shocking 
induced mass loss is not self-limiting, but continues at a con- 
stant (orbit-averaged) mass-loss rate without limiting the cluster 
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to its (perigalactic) tidal radius. This behaviour, which is at odds 
with previous results by Gnedin et al. (1999), is probably a con- 
sequence of the eccentric orbit, which enables the cluster to set- 
tle into a new equilibrium near apogalacticon, i.e. unaffected by 
the perigalactic tidal field. 

Our findings imply that initially not tidally limited clusters 
of low to medium concentration and moving on eccentric orbits 
may never become tidally limited and, hence, can be quickly de- 
stroyed even when their relaxation time is long. Clearly, more 
theoretical and numerical work is needed to clarify and quantify 
the issues of cluster destruction by tides, in particular the effect 
of the Galactic orbit for the relative importance of disk shocks 
in the evolution of globular clusters. 

Palomar5 is an excellent example for a super-tidal globular 
cluster which extends to four times its perigalactic tidal radius. 
Some of our simulations actually reproduce this property to- 
gether with the observed very low velocity dispersion. Both 
are the immediate consequences of the last disk shock, which 
heated the cluster and made it expand to its present size. Our 
simulations also uniquely predict the destruction of Pal 5 at its 
next disk shock in ~ 110 Myr. After that the cluster will loose 



its identity and eventually only an increasingly separated pair of 
tidal tails will remain. 

The tidal tails observed within a few degrees from Pal 5 
(paper II) show substantial structure which cannot be repro- 
duced by our simulations. We argued that this structure must 
be due to Galactic sub-structure (not accounted for in our sim- 
ulations), such as giant molecular clouds, spiral arms, or even 
dark-matter sub-halos or massive compact halo objects. More 
work is needed to quantitatively address the origin of this struc- 
ture. 

The tidal tails may also be used to infer the mass of the Milky 
Way within 18.5 kpc to an accuracy of <~ 10% by measuring the 
radial velocity or distance modulus of the tail stars over a part of 
the tidal tails as large as possible. To this end, an extensive ob- 
servational programme is needed to uncover the further extent 
of the tidal tail beyond the SDSS data. 



We thank the anonymous referee for his quick yet useful and 
detailled report. EKG gratefully acknowledges partial support 
from the Swiss National Science Foundation. 
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